Entanglement at the quantum phase transition in a harmonic lattice 
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The entanglement properties of the phase transition in a two dimensional harmonic lattice, similar 
to the one observed in recent ion trap experiments, are discussed both, for finite number of particles 
and thermodynamical limit. We show that for the ground state at the critical value of the trap- 
ping potential two entanglement measures change abruptly. Entanglement thus indicates quantum 
phase transitions in general; not only in the finite dimensional case considered in Q. Finally, we 
consider the thermal state and compare its exact entanglement with the result of the temperature 
entanglement witness, introduced in 0). 



INTRODUCTION 

Coupled harmonic chains with short and long range 
interactions are ubiquitous in science and engineering. 
Their application to calculate the phononic heat capac- 
ity by Einstein 0] marks the birth of solid state physics. 
Beyond physics, harmonic chains feature in chemistry 
and biology, where they are used to model behaviour of 
macro-molecules, such as DNA jj| and cell membranes 
In the last decade harmonic systems have been re- 
vised using techniques developed in quantum information 
science to study correlation properties in the quantum 
regime and particularly at small temperatures 0, 0, B Q • 
Thermodynamics has been very successful in characteris- 
ing "standard" phase transitions that occur at finite tem- 
perature when a macroscopic parameter, such as pres- 
sure, is changed [loj ]. Quantum phase transitions (QPTs) 
appear at zero temperature ll| and are due to the change 
of an external parameter, such as the trapping potential 
of an ion trap. These transitions are driven by quantum 
fluctuations and have been linked to entanglement for the 
case of finite dimensional systems [H [HJ ■ 

In this paper we study a QPT in a continuous variable 



system [13|, ll4| - a system of trapped ions interacting via 



a long-range Coulomb repulsion, see Fig. [TJ The ions are 
trapped by two external potentials and align in a linear 
configuration for big enough transversal trapping poten- 
tial, v t . When v t is decreased the system undergoes a 
phase transition and the new equilibrium state forms a 
zig-zag configuration. This model is motivated by ion 
trap experiments 15|j li|, 17 1, where such a QPT occurs 



[ 19j]. Recent analytical studies of the transition us- 
ing Landau theory |20| allowed to determine the system's 
classical behaviour at the transition point. Moreover, the 
numerical treatment of the quantised system of a few ions 
promised the possibility of simulating linear and nonlin- 



ear Klein-Gordon fields on a lattice [2l|. However, a com- 
prehensive analytical study of the quantised system has 
so far been lacking due to the complexity of the system. 



Here we model the ion trap scenario as a lattice of har- 
monically coupled oscillators and present a first quanti- 
tative characterization of the entanglement inherent in 
both ground state configurations of the ions. For fi- 
nite dimensional systems QPTs of first (second) order 
are characterised by a discontinuity (a discontinuity in 
or divergence of the first derivative) of the negativity [lj . 
We show that also in this continuous variable system the 
structure of the entanglement, measured by the negativ- 
ity and the von Neumann entropy, changes abruptly at 
the critical point and indicates the occurrence of a QPT, 
in a similar way as classical correlations indicate standard 
phase transitions. The first derivative of the negativity 
between two neighbouring ions has a finite discontinu- 
ity and the von Neumann entropy of a single ion has a 
divergence in the first derivative. 



The quantum fluctuations that cause the QPT are 
most dominant at zero temperature. However, as exper- 
iments are performed at small, but finite temperatures 
it is important to know how temperature affects these 
fluctuations 



22j. In the final part of our paper we dis- 
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cuss thermal states and find that the sharpness of the 
QPT, indicated by the entanglement, fades out with in- 
creasing temperature. Another macroscopic consequence 
of quantum fluctuations is the lowering of the energy of 
the system We compute up to which temperature 
the thermal state has a lower energy than any separable 
state. This temperature witness could be implemented in 
an ion trap experiment by measuring the average energy 
(i.e. mean excitation). Although the model is motivated 
by ion traps realisations include the vibrational motion 
of molecules (nuclei in the electronic potential) and 
optical lattices [23(. 
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v t > v t . 



where V 3 ■ = \ 
of the sites with ■ 



Q 



k^j \fk-fj\ 



are the Coulomb potentials 



VjiVj) absolute coordinates of the 
sites and Xj t o are the equilibrium positions in x direction. 
N is the number of particles, Q is the charge of the ions 
and v (y t ) are the trapping potentials in x (y) direction. 



We assume periodic boundary conditions, 



Tj = r j+N . 




Vt < Vtfivit 



FIG. 1: Sketch of the harmonic lattice under consideration: 
Each site is trapped by two external potentials, v in x di- 
rection and vt in y direction. For clearness only the nearest 
neighbour coupling is indicated, however all ions interact via 
a long-range potential which we approximate harmonically. 
This could be for instance a Coulomb potential such as is 
common in ion experiments. The equilibrium distances be- 
tween the sites are assumed to be equidistant, with lattice 
constant a in ^-direction and b along the ^/-direction. In /, 
the transverse trapping potential vt is larger than the critical 
value and the ions arrange linearly. As displayed in II, de- 
creasing the transverse trapping potential below the critical 
value vt,crit leads to a QPT causing the ions to move outwards 
and form a two-dimensional zig-zag structure. 



To calculate the entanglement measures we approxi- 
mate the Coulomb potential to second order and expand 
about the equilibrium positions. The key step is then 
to diagonalise the Hamiltonian into a set of uncoupled 
modes, the lattice vibrations, with which analytic expres- 
sions for the measures can be obtained. Similar to the 
classical calculation jipj . we use a simplified model with 
equidistant equilibrium position in x direction, spaced by 
the lattice constant a. Such condition can be realised for 
the central ions of a long ion chain inside a linear Paul 
trap 2m or for ions confined in a ring of large radius 



25] 



For big trapping potential the sites are arranged on a 
single line, i.e. Xj t o = aj and j/j.o = 0, while for small 
enough v t , the equilibrium positions become Xjo = aj 
and ijjfi = (— 1) 2 and a two-dimensional zig-zag configu- 
ration is formed. The equation determining the deviation 
b in y direction is obtained by summing the linear terms 
over all sites and requiring it to vanish, 



THE MODEL 



We consider the Hamiltonian, see Fig. [TJ 

N / „2 1 „2 2 2 

h m = E I 2m + — - + -r-yj + V i 



3=1 



\m V 2 t =Q 2 J2 



T=2 /+i Vr 2 a 2 + b 2 ' 



(2) 



where r = k — j numbers the neighbours of each sites. 



(1) The Hamiltonian becomes 
I 



N / 2 , 

= / p XJ +p: 
^ ' 2m 



mv. 



t „,2 



Q 2 N 



— ( rf r ( x j ~ Xj+rf + d v T {yj - y J+T ) 2 + d x T y (xj - x j+T ){ yj - y j+r )) , 



J=1,T>0 



with Xj = Xj ~ stj.o and yj = jjj — y^o the deviations from 
equilibrium. Furthermore, the d x ^ yjXy denote the second 
order Taylor coefficients of the Coulomb potential which 



are, for the linear and zig-zag configuration, 



dl 



1 



(at) 3 
1 



2(ar) 3 



and dl 



2r 2 a 2 - S T:0dd b 2 
V(™) 2 +$r,oddb 2!t 



d% v = Q and d% v = 5 T a , 



2yJ{Ta) 2 + 5 TMd b 2 " 
{-lySrab 



dd ■ 



(3) 



25 T ,oddb 2 - T 2 a 2 
and dl = ^ (4) 



2yJ(ra) 2 +b r 



.(5) 
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CALCULATION OF ENTANGLEMENT 
MEASURES 

We are interested in the behaviour of the entanglement 
between the sites in the chain for varying transverse trap- 
ping potential, v t , and particularly at the point of crit- 
icality, ft, C rit- We calculate the negativity, En, for two 
modes regardless of all others, for instance, the x degrees 
of freedom of two (neighbouring) sites. To measure the 
correlation of one mode with all other modes, for instance 
the entanglement between the y degree of freedom of a 
single site with all the other degrees of freedom in the 
chain, the von Neumann entropy, Sv, is used. The same 
measures could be used to evaluate the entanglement for 
other splits of the available degrees of freedom, see @ 
for block entanglement. Other entanglement measures 
are available, such as the entanglement of formation [2(| . 
Yet they are very hard to calculate in this continuous 
variable scenario and we will be content with the two 
measures as stated. To see the effect of the long-range 
interaction we compare the full long-range (LR) Coulomb 
Hamiltonian with a cut-off version in which only nearest 
neighbours interact (NN), similar to 0. 

For the linear configuration the Hamiltonian decouples 
into x and y part. A discrete Fourier transformation for 
the x (similar in y direction) of the form 



1 



N 



Xi 



Pxj 



1 N 



xl 



(6) 



(7) 



maps the space coordinates of the sites into diagonal 
modes, the lattice vibrations or phonons. The diag- 
onal Hamiltonian is then H = Y^iLi ^xi (fi x i + 5) + 



where 



n x i 



p^Pxi mu> xl x; X, 1 

2hmto x i + 2h 2' 



J2i=i ^yi [ n yl ' 2 J 

and similarly n y i, are the number operators in x and y 
direction for mode I, and the frequencies are 



sin 



' Tdr_ \ 

v N J 



T>0 



sin 



' ttIt \ 
. N ) 



(8) 



(9) 



T>0 



where C = -^r. The asymmetry of the system is rc- 
fleeted in the dispersion relations. ui x i is always real, 
whereas ui y i would become complex for small values of 
transverse trapping potential vt . This is where the quan- 
tum phase transition occurs. The critical value vt,crit for 
LR interaction is Vt,crit ~ v0l3C while for NN interac- 
tion it is Vt.crit = V0.5C. 

To diagonalise the Hamiltonian in the emerging zig-zag 
configuration, where the previously independent phonons 



now couple in x — y-dircction, see Eq. ([3]), we need to 
amend the transformation in the y direction Eq. ([6]) and 
transform with 



y = ' y e i%j(,l+N/2) Y 

3 Vn Y 



Pyj 



l —y e -^ l+N ^p yl . 
'iv y 



(10) 

(ii) 



The additional factor e™- 7 = (— l) 3 compensates the al- 
ternating sign of the x — y coupling in the Hamiltonian. 
Expressed with a coupling matrix 



Mi 



m~2 



m ~ 



-US 





1 

2m 

rn~2 



xyl \ 







V 







2^1 







(12) 



J— 

2m / 



the Hamiltonian can now be written as 



N 



with 



1=1 



coefficients 



[ X t \ 

Yi 

\PyJ 



(13) 



x, y and 



■>xyl 



= 91 v 



3rab 



^sm^lr/N). 



>0 < odd 2y/(ray+^° 

For each I we need the symplectic transformation [27| , 



denoted by Si, that diagonalises the matrix Mi, 



Si Mi SI = diag 



2 ' 2 



(14) 



where 



+ tfi sl^li - <^) 2 + 4 C ( 16 ) 



are the symplectic eigenvalues of Mi. These transforma- 
tions are given by Si = 



( 



















I = ^xl 



2u> X yly 








(17) 



+ yj ("Ii ~ tfi) 2 + ^lyi and <Pi = 



The new normal modes are 



p 

Wl 



= S7 



/ Xi \ 

Px, 

Yi 

\ p yJ 



(18) 
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2h 



and using the number operators n v> i 

— i, and similarly for w, we find the fully diago- 
nalised Hamiltonian for the zig-zag configuration H = 

J2i [^>v,t (n v ,i + |) + h^w.i (n w ,l + 5)] • 

For calculating the negativity we use the PPT criterion 
for bipartite Gaussian states [28|. Any bipartite state 
whose partially transposed covariance matrix does not 
fulfill the uncertainty relation is entangled. For two sites 
in two dimensions the criterion requires the symplcctic 
eigenvalues of a 8 <g> 8 covariance matrix. However, all 
expectation values coupling the x and y direction, i.e. 
(xj yj), {p X j p y j) etc., vanish so that two sets of each two 
conditions, as used in emerge. For each site j these 
separability conditions are 

< S x>2 {v tl r) = ± ((x ±x 0+T f) ((p Xj T^ +T ) 2 ) ~ 1 (19) 

and similarly for the y direction. The expectation val- 
ues needed here can be calculated using the transforma- 
tion rules into the diagonal modes, Eq. ((6]), Eq. (fT0|) and 
Eq. (fT8|) . If one of the inequalities is violated then entan- 
glement exists between the j-th site and its r's neighbour 
and the negativity, 



and replace the sums with integrals 

1 r /2 

^ J-tt/2 



£1,2(^,7") = j u> x , a ( 1 ± cos ( z ^- r r ) ) da 



2na 
N~ 



r /2 (l T coster)) 

/ ^ - 1 

J-n/2 w x,/3 



(22) 



These integrals can be evaluated numerically and the 
plots are shown in Fig. [2] (red-solid line). At the criti- 
cal point vt,crit = % for NN interaction the negativity 
can be calculated analytically as 

4 r /2 fc 



cos(a) (1 — cos (2ar)) da 



/* (1 + cosW)^ 





4 2 
^2 ■ 3 



e ■ cos(/3) 



2 - 1 w -0.46 



(23) 



which leads to a value of the negativity of En ~ 0.308. 
In a similar fashion, the symplcctic eigenvalue r and the 
von Neumann entropy for both, x and y direction, can 
be evaluated in the thermodynamical limit. 



En = ^2 max 0, - In S k + 1 



(20) 



fe=i 



measures their degree of entanglement. The two criteria 
Si and S*2 witness two types of entanglement. For ex- 
ample, the EPR pair originally considered in [2^] shows 
violation for 5*2 but not S\. 

The von Neumann entropy of a single site j in either x 
or y dimension is obtained following [301 ] with the formula 



S v (rj) 



^Inf^V^m^- 1 



(21) 



where rj = (x'j)(p'^ ~), and similiarly for the y direc- 
tion, is the symplectic eigenvalue of its reduced state. 

Thermodynamical limit (N — » 00) 

To obtain the negativity for the x entanglement at zero 
temperature, we evaluate Eq. p9[) in the ground state 
leading to 



i.k 



and similarly for the entanglement in y direction. In the 
thermodynamical limit we substitute ^ = a, = /?, 



BEHAVIOUR OF ENTANGLEMENT AT ZERO 
TEMPERATURE 

Fig. [2] displays both entanglement measures for de- 
creasing transverse trapping potential v t . In the upper 
plots (a and b) both measures for the ^-entanglement 
are constant in the linear regime. This is because the 
phonons in x direction are independent of the trapping 
in y direction. At the critical point both negativity and 
entropy are not differcntiablc. Decreasing the trapping 
potential beyond the critical value, where the zig-zag con- 
figuration is formed, the rr-cntanglcmcnt is reduced due 
to the emerging x — y coupling. Because each site in (b) 
couples to several different neighbours, the entanglement 
between two distinct sites disappears. 

In the lower plots (c and d) both entanglement mea- 
sures for the y-entanglcment grow with decreasing trap- 
ping potential v t in the linear configuration. The even 
numbered ions oscillate exactly out of phase with the 
odd numbered ions due to the repulsive Coulomb poten- 
tial. The smaller is t , the larger these quantum fluctua- 
tions around equilibrium position become. At the critical 
point the fluctuations become strong enough for causing 
the ions to move outward. The entropy diverges and 
the negativity reaches its maximal value of En f=a 0.308 
where it is not differentiablc. 

For the nearest neighbour coupling (a and c) the nega- 
tivity and entropy show qualitatively the same behaviour. 
This can be understood easily as any entanglement of 
an site with the rest of the chain, measured by the von 
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FIG. 2: This graphic shows the von Neumann entropy of a single site and the negativity between two neighbouring sites both, 
for x-entanglement (upper plots) and y-entanglement (lower plots) of the ground state (T = 0) . We compare the case where 
only nearest neighbour interact (left plots) with the full long-range Coulomb Hamiltonian (right plots). The numerical values 
for the plots are for the ground state (T = OK) in the thermodynamical limit (N — > oo) and Q — 1, ra = 2. For NN the lattice 
constant is set to a = 1, and for LR it is a = 14/15 and the plots include interactions up to the forth neighbour. The change 
of entanglement at the critical point, indicated by the vertical line, is clearly visible in all four plots. 



Neumann entropy, is created by the coupling with only 
the nearest neighbours. For the long-range interaction (b 
and d), where a single site couples to all other sites, there 
are significant differences between the two entanglement 
measures. While the negativity of two nearest neighbour 
sites vanishes after a threshold value of v t , each single 
site remains entangled with the rest of the chain, as seen 
by the positive value of the von Neumann entropy. 

Both entanglement measures are functions of the eigen- 
values of the Hamiltonian. Therefore abrupt changes in 
entanglement can signal non-analyticity of ground state 
energy, which is associated with QPT's. In [lj it was 
shown for finite dimensional systems that (under certain 
conditions) a discontinuity in or divergence of the first 
derivative of the negativity is both necessary and suffi- 
cient to signal a QPT. It seems intuitive that a similar 
characterisation holds also for continuous variable sys- 
tems. Here the divergence of entropy and finite discon- 
tinuity of the first derivative of negativity shows a QPT 
of second order. After the phase transition the violated 
negativity criterion switches from S\ to 1S2. Lowering 



the trapping potential further leads to another "critical" 
point c y at which the negativity becomes zero while the 
entropy reaches its minimal (non-zero) value. 

The additional critical point c y (and c x in (b)) only ap- 
pears when the interaction is harmonic as is the case in 
our second order approximation of the Coulomb model. 
The point is due to a sign change of the second order 
coefficient d% {d%) at c y (c x ). As a consequence the in- 
teraction switches from repulsive to attractive. When 
anharmonic terms are taken into account, as in the nu- 
merical treatment in plj . these points vanish. An intu- 
itive way of understanding the switch in the harmonic 
approximation is as follows, see Fig. [3l Decreasing the 
trapping potential changes the relative strengths of the 
inner and outer potential for the motion in the y direc- 
tion. This leads to the change in phase in the relative 
momenta of neighbouring sites. In one configuration the 
relative potential favours momenta in the opposite direc- 
tions, while the other configuration favours motion in the 
same direction. This is reflected in the change from Si to 
S2. However, although the negativity vanishes at these 
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WITNESSING ENTANGLEMENT AT FINITE 
TEMPERATURE 

One consequence of the entanglement of the sites is a 
lowering of the energy of the system 0] ■ This can be seen 
by assuming that the thermal state of the system 



N 



PS = 



y tr[ e -^]' 



(24) 



FIG. 4: This graphic shows the negativity in y direction of 
two neighbouring ions for decreasing trapping potential and 

increasing temperature (in units of [u t ] = <J and [T] = 



h 
2k H 



N = 20, 



a = 1, m = 2). The red line indicates 
the critical temperature, obtained with the energy witness 
argument, below which entanglement of some form must be 
present in the chain. 



with p 



k B T 



the inverse temperature, is separable, 



i.e. decoupled between modes. Then an effective, sin- 
gle site Hamiltonian can be obtained by removing all 
second order couplings between the different sites, i.e. 



(( x j 



.r 



3+tJ 



IPs 



etc., and then H c s 

„2 



3 ' J+T 

Y] ■ Hj with 



\-r)ps — "\-°j/ps 



with 
- - j 



V 



yj 



2m 



+ f (n^J + Sfe 2 + O^ClVl) (25) 



\ 



— dS. and Q v = 

T=£0 

and Cl~,. - 



\ 



T#0 



— E d Z (26) 

TO ' J 

T=iO 



Using the transformations Eq. ©, Eq. (|TD|) and Eq. (fT5)) . 
(H e g) can be fully diagonalised. It then becomes clear 
that the internal energy U = (i? ft) f° r any separable 
state is bounded from below by 

Nh 

(H eS ) se p >—{ttx+ My + £l X y) (27) 

and any state having a smaller energy must be entangled 
between the individual sites. 



We now want to see how the ground state situation 
is modified at non-zero temperature. As the energy of 
the thermal state, i.e. the mean excitation of phonons, 
increases with temperature there exists a critical tem- 
perature, T c , for each value of the trapping potential at 
which the thermal state matches the energy bound. The 
negativity between the y degrees of freedom of two near- 
est neighbours for NN Hamiltonian is evaluated numer- 
ically and plotted in Fig. [3] where the critical tempera- 
ture for full separability is also indicated as a red-line. 
When the trapping potential is lowered, the negativity 
increases until reaching its maximal value at the criti- 
cal vt,crit- Further lowering leads to a decrease of the 
negativity until it vanishes at point c y . However, when 
further decreasing v t the negativity grows again, yet the 
two entanglement criteria S\$, are now switched. As ex- 
pected, increasing the temperature leads to smaller val- 
ues of negativity and smoothens the entanglement mea- 
sures to make it diffcrcntiablc at the critical point. For 
large v t the negativity is small, but remains finite un- 
til relatively high temperature. The sharp peak at the 
QPT remains almost constant for finite T and decreases 
fast. Thermal states within the red outlined area have 
a smaller energy than any separable state and their en- 
tanglement is therefore detected by the energy witness. 
Remarkably, states at c y are entangled for temperatures 
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up to T c [l/v t = c y ] = 0.12^^, even though there is no 
nearest neighbour entanglement in y direction. This is 
because there is still nearest neighbour entanglement in 
x direction and possibly also multipartite entanglement 
the chain. 



CONCLUSIONS 

In this paper, we revised the classical phase transi- 
tion in a long range harmonic chain [20| using a fully 
quantised model. Two measures of entanglement display 
critical behaviour: The von Neumann entropy diverges 
at the critical point while the negativity is not differ- 
cntiablc. Thus also in this continuous variable system 
entanglement indicates a QPT, as previously shown for 
discrete systems Furthermore, our results confirm 
that this phase transition is of second order as indicated 
in [20l | . At finite temperature, the negativity still dis- 
plays critical behaviour, as seen in Fig. |4j however the 
non-differentiable cusp fades out quickly with increasing 
thermal noise. Tuning across a QPT provides a means of 
changing the amount and structure of continuous vari- 
able. Experiments with ion-traps are ideally suited to 
study QPTs with great precision. 
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